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1.0 INTRODUCTION 

Previous error analyses conducted by the Guidance and Dynamics 
Branch (GDB) of NASA have used the Guidance Analysis Program (GAP) 
as the trajectory simulation tool. Current plans are to conduct 
all future error analyses using the Space Vehicle Dynamics Simula- 
tion (SVDS) program. A study has been conducted to compare the 
Inertial Measurement Unit (IMU) error simulations of the two pro- 
grams. This paper presents results of the GAP/SVDS comparison and 
defines problem areas encountered while attempting to simulate IMU 
errors, vehicle performance uncertainties and environmental un- 
certainties using SVDS. An evaluation of the SVDS Linear Erro»" 
Analysis (LEA) capability is also included. 

2.0 DISCUSSION 

2 . 1 GAP/SVDS Comparison for IMU Error Sources 
To evaluate an IMU simulation, the effects of guidance/navigation 
interfacing should be considered. The Mathematical Physics Branch 
(MPB), Software Development Branch (SDB) and GDB previously compared 
SVDS and Navigation Analysis Program (NAP) simulations for IMU 
errors. (See Reference 1.) However, the NAP and SVDS comparison 
was made at trajectory times prior to Solid Rocket Booster (SRB) 
staging. The effects of SVDS program phasing and guidance/ 
navigation interfacing v/ere not considered. This GAP/SVDS comparison 
considers both of these effects by making comparisons of trajectory 
data at main engine cutoff (MECO). 
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2.1.1 Selection of Nominal Trajectory and IMU Error Source 
Magnitudes 

The last error analysis conducted by GDB was in June 1974. The 
trajectory used as a basis for the IMU error analysis vas a Base- 
line Reference Mission 3A Abort Once Around (AOA) boost profile. In 
order to obtain comparable trajectory data from the IMU error 
sources, an SVDS simulation was developed to match this GAP nominal 
trajectory. A list of some key trajectory parameters of the GAP 
nominal trajectory and the SVDS trajectory is given in Table I. 

The table contains a comparison of the GAP and SVDS conditions at 
main engine cutoff (MECO). The trajectory differences at MECO are 
insignificant indicating that the GAP and SVDS trajectories compare 
closely enough for use in this analysis. The indicated SVDS 
trajectory is used as the basis for the SVDS IMU error analysis. 

IMU error sources as defined in the last GDB error analysis were 
selected for use in this study. These values may not be the most 
recent evaluation of the uncertainties but were selected to obtain 
valid comparison data. Table II contains the error sources and 
their 3-sigma uncertainties in both GAP c.id SVbS input units. The 
SVDS simulations use the 3-sigma uncertainties for the error 
sources. 

2.1.2 Trajectory Data Comparison 

Table III contains the GAP/SVDS trajectory comparison data. Data 
are presented for deviations in position and velocity at MECO due 
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to each IMU error source simulation. The deviations are a result 
of 3-sigma uncertaintines and are computed as 

L - (actual perturbed state vector component at MECO) 

- (nominal state vector component at MECO) 

and are in an Earth-Centered Inertial (ECI) coordinate system. For 
each error source, deviations from the nominal are presented for 
both GAP and SVDS simulations. GAP data is the first line of state 
vector deviations for each error source and the SVDS deviations 
are presented in parentheses. 

Examination of the trajectory deviation data reveals the following 
differences in GAP and SVDS simulations: 

a. For some error sources, absolute value of the deviation in each 
state vector component is approximately the same but there is 

a difference in the signs of the deviations. 

b. For some error sources, a difference exists in the definition 
of spin axis and output axis components. 

c. For some error sources, the comparison between GAP and SVDS 
deviations show a large percentage variation in Z component of 
velocity. 

An effort was undertaken to investigate these differences. 
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2. 1.2.1 Sign Differences in GAP and SVDS Deviations 
Examination of the programming code for GAP and SVDS reveals that 
sign differences exist in how the programs model some of the IMU 
error sources. The modeling differences are reflected by corres- 
ponding trajectory deviations having similar magnitudes but opposite 
signs. The modeling differences result from the fact that the IMU 
model simulated by GAP is defined according to the conventions 
established for the Apollo project while the SVDS IMU model is 
consistent with current IMU definitions established by MPB (Reference 1). 
The following error sources are modeled with different sign conventions 
in GAP and SVDS and the differences are indicated by the deviation data 
in Table III: 


a. free gyro drift bias (Z IMU error) 

b. gyro spin axis acceleration sensitive drift (X and Y IMU errors) 

c. gyro output axis acceleration sensitive drift (Z IMU error) 

d. accelerometer input axis misalignment toward the output axis 
(X and Z IMU errors) 

e. accelerometer input axis misalignment toward the spin axis 
(Y component). 


2. 1.2. 2 Axis Definition Differences 

An additional difference between GAP and SVDS exist in their 
definitions of spin axis and output axis for some of the error 
sources. The program definitions of spin axis and output axis are re 
versed when considering some components of IMU accelerometer input 
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axis misalignment and gyro acceleration sensitive drift. In 
determining how sensed velocity is to be perturbed for the Y 
component of the IMU error source (accelerometer input axis mis- 
alignment) » GAP perturbes the Y component of sensed velocity by 

a. Vy = Vy + (error source magnitude) 
for misalignment toward spin axis and by 

b. Vy = Vy + V^* (error source magnitude) 

for misalignment toward output axis. The SVDS simulation models 
the spin axis and output axis oppositely. Equation (a) is used 
for perturbed sensed velocity for output axis and equation (b) is 
for spin axis in SVDS. GDB concurs with the MPB (SVDS) convention 
for these error sources. The GAP/SVOS comparison data of Table III 
is arranged to the GAP definitions for uniformity of comparison. 

Similarly, for gyro acceleration sensitive drift (Z comporent of 
the IMU error), GAP models the perturbed sensed velocity X and Y 
components as 

2 

®1’ ^ '^Y * source magnitude) 

~ ^Y ^ '^X ^Y* source magnitude) 


for spin axis and as 


h* ^ '^X '^Y* source magnitude) 

2 

bg* Vy * Vy + Vy * (error source magnitude) 

for output axis. The SVDS simulation models the spin axis and out- 
put axis oppositely. Equations (a-|) and (a 2 ) are used for output 
axis and equations (b-j) and (b 2 ) are used for spin axis when con- 
sidering the Z component of IMU error in gyro acceleration sensitive 
drift. Table III is arranged according to the GAP definitions for 
this error source for uniformity of comparison. 

2. 1.2. 3 Z Velocity Deviation Differences 

After making allowances for the previously discussed sign and axis 
definition differences, examination of the GAP and SVDS-generated 
state vector deviations of Table III shows that the only state 
vector component for which large percentage variations exist between 
GAP and SVDS is in the Z component of velocity for some of the 
error sources (e.g. gyro drift bias or accelerometer bias). 

MECO conditions were investigated to see if the variation is a 
result of a variable MECO time or the inaccuracy of the cutoff 
velocity magnitudes. Neither of these proved to be the case. 

All of the problem cases had accurate cutoff times and comparable 
MECO velocities. 

To determine which set of Z velocity deviations is more plausible 
(GAP or SVDS), an attempt was made to correlate velocity 
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deviations at MECO and position deviations. Consider the state 
vector deviations given in Table III for gyro input axis acceleration 
sensitive drift. Assuming that the X and Y components of velocity 
deviation are constant from liftoff to MECO in the sample case, it may be 
observed that a 2.5 ft/sec velocity deviation in X projects 

into a 460 ft deviation in position. Similarly, a -.04 ft/sec Y 
velocity deviation results in a -7 ft deviation at MECO. (See 
Table III.) Determining a ratio of position error to velocity 
error yields 

X position _ Y position ^ 180 ft 

X velocity " Y velocity ft/sec 

for the indicated error source. Now consider the Z velocity devia- 
tion. Using the GAP value and the indicated position to velocity 
ratio yields 

Z position . 180 * -2 = -36 ft. 

Using SVDS values yields 

Z position . 180 ♦ -7 = -126 ft. 

The actual position deviations for the error source are -127 feet in 
GAP and -123 feet in SVDS. The Z velocity and resulting position 
deviations show better correlation to the X and Y components in the 
SVDS. Similar results may be obtained for other sources- 

The programming code 'of GAP and SVDS were examined to resolve the 
Z velocity deviation differences. No errors were found in either 


- 7 - 


program, The lack of comparison between GAP and SVDS for this 
velocity component remains an unresolved problem. However, since 
the magnitude of the Z component of velocity deviations are small 
and since the SVDS values appear to be acceptable, no further 
Investigation of this discrepancy was attempted. 

2.1.3 Conclusions 

State vector deviations resulting from SVDS simulation of IMU errors 
are comparable to the deviations resulting from GAP simulations. 

SVDS may now be used as a simulation tool for generating an IMU 
error analysis. The sign and axis differences existing between GAP 
and SVDS present no problem since MPB has confirmed the IMU model 
within SVDS. (See Reference 1.) 

2.2 Checkout of SVDS Dispersion Analysis Capability 
Besides simulating IMU errors, SVDS must be able to simulate vehicle 
performance, aerodynamic and environmental uncertainties if it is 
to be an effective dispersion analysis tool. To checkout SVDS 
capabilities, trajectory simulations were developed using 3-sigma 
uncertainties in each of the following; 

a. vehicle vacuum thrust, specific impulse and propellant loading 
for both the SRB and main engines. 

b. inert weight for the SRB and external tank and orbiter 


c. axial force coefficient and base drag 


d. Vandenberg hot and cold atmospheres. 


The resulting MECO insertion weight deviations were comparable to 
previous GDB dispersion analysis results. A 6AP/SVDS comparison 
of trajectory and weight dispersions is not made in this report 
because: 

a. since the last GAP error analysis was conducted, new SRB perturba 
tion techniques have been defined (Reference 2). 

b. nominal SVDS programming logic has previously been checked. 

2.2.1 Problems Encountered in SVDS Dispersion Analysis Simulations 
The following SVDS problems are encountered when attempting dis- 
persion simulations: 

a. stacked cases (multiple trajectory simulations in one job 
submittal) result in overflow of allotted tabular input 
capability. 

b. Vandenberg hot and cold atmosphere models do not execute. 

Simulation of SRB thrust and specific impulse uncertainties require 
input of several sets of tabular data (SVDS input parameter TABLE). 
When attempting stacked cases for these siri’ul ations , overflow of 
the TABLE length is encountered. According to SVDS documentation, 
input tabular data is overlayed for stacked cases. The over flow 
of TAL.LE encountered for tfiese cases indicates ttiat in multiple 
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SVOS cases, tabular data Is added In series not overlayed. Soft- 
ware Development Br;«.ich (SDB) was notified of the problem and 
Indicates that the program alterations required to fix this problem 
would be Included In the next SVDS update; however, they do not 
plan to attempt any checkout cases to verify the program modifications. 

It should be noted that one program fix for this problem Is to In- 
crease the dimension of allowable tabular Inputs. This was not 
attempted In this study because of a lack of unused core on the cur- 
rent version of SVDS. 

Attempts were made to simulate the Vsneenberg hot and cold atmospheres 
as environmental uncertainties. Neither of these SVOS options will 
execute. The programming code exist for modeling the uncertainties 
but program modifications are necessary before SVDS will execute the 
options. A discrepancy report (Reference 3) has been submitted to 
SDB indicating this problem. 

2 . 2.2 Conclusions of SVDS Mspersion Analysis Checkout 
To efficiently conduct a dispersion analysis, the trajectory 
simulations for t:.e uncertainties should be generated in a limited 
number of computer job submittals. SVDS cannot currently handle 
the required stacked cases which involve tabular data. Atmospheric 
uncertainties are also not executable on SVOS. SDB support is 
required in both these areas before a complete dispersion analysis 
can be conducted. 
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2.3 Checkout of Linear Error Analy s is (LEA) Capability 


The current directive from GDB Indlcdtes that dispersion analyses 
will be conducted as single error source cases. A LEA technique 
win be used to statistically correlate trajectory deviations for 
all uncertainties by developing covariance matrices at various 
flight events. As part of this study, an effort was undertaken 
to exercise and verify the LEA capability In conjunction with SVDS. 
The study was hampered by the lack of good documentation defining 
techniques for use of the LEA processor. 

2.3.1 LEA Procedures 

The LEA option of SVDS determines the following output data: 

a. trajectory parameter deviations defined In a local horizontal 
coordinate system (LHS) 

b. the root- sum- square (RSS) of the trajectory parameter devia- 
tions 

c. a covariance matrix relating all of the error sources. 

To exercise the LF.A processor, SVDS trajectories must first be 
simulated. For the LEA processor to generate covariance data and 
RSS data for the desired trajectory parameters at the desired 
flight events, the fo1 leaving procedure Is required In developing 
the SVDS trajectories: 

a. determlno the state vector conponorits and performance parameters 
for which I TA results are desired 
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b. determine the flight events at which LEA results are desired 

c. generate a nominal SVOS trajectory and trajectorua for each 
of the error sources using the 3-s1gma uncertainty 

d. c ate a tape containing tie perturbed trajectory parameters 
at each of the required flight events for all the SVOS 
trajectories. 

The SVOS- generated tape should contain a file for each SVDS trajectory 
(a nominal followed by the perturbed trajectories). Each file of the 
tape should contain records for only the flight events required for 
the LCA. Each record should contain the trajectory parameters to be 
used in the LEA. Schematically the tape may be represented as; 

File 1 (Nominal) 

Record 1 Trajectory parameter data 

Record 2 Trajectory parameter data 


t.d of file 

File 2 (Perturbed Case) 

Record 1 Trajectory parenicter data 

Record 2 Trajectory parar;oter data 


Fnd of file 

File i (PerLurlKjd Case) 
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To generate a SVDS tape written as depicted, the following SVDS 
Inputs should be used: 

a. PLOT = The desired trajectory parameter data 

b. IPLOT = 10000 In those SVDS program phases which terminate 

at a flight event specified for LEA 
IPLOT =0 In all other SVDS program phases 

c. MCPT = 1 

It should be noted that the SVDS tape write and LEA tape read units 
cu.rently do not match if default values are used in both 
processors. SVDS assigns the tape write to unit F while the LEA 
reads unit D. This poses no problem if the SVDS trajectory simula- 
tions and LEA processor are submitted separately. 

The tape assign control cards (ASG = unit i) must be identified as 

unit F (in SVDS) and unit D (for LEA processor). If the SVDS 
trajectory simulations and LEA processor are run in one job sub- 
mittal, the output and input units must be matched. This can be 
controlled by inputting IPUNIT = 4 in SVDS (forces output onto 
unit D) £T setting JUNIT = 8 in LEA inputs (forces reading of unit 
F). The tape assign cards should be input accordingly. 

As previously indicated, the LEA processor uses trajectory disper- 
sions developed using 3-sigma uncertainties in the error sources. 
The LEA processor reads the SVDS state vectors of the nominal 
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trajectory (developed in an earth-centered inertial coordinate 
system) and forms a LHS based on the nominal data. The LHS 
coordinate system is defined by: 

0 = iT/m 

V = (R X V X R)/lR X V X R! 

M = U X V 

where ^ and Y are the radius and velocity vectors of the nominal 
trajectory at the trajectory time for which the LEA exercise is 
specified. The LEA processor then rotates the state vector devia- 
tio. ? of the perturbed cases into the LHS and combines the devia- 
tions by a root-sum-square (RSS) process. The resulting RSS is 
based on 3-sigma uncertainties. The LEA processor then forms a 
covariance matrix indicative of all the simulated error sources. 

The covariance matrix is based on 1-sigma level of confidence for 
the error sources. 

2.3.2 LEA Difficulties Encountered 

Other than the lack of documentation defining LEA procedures, the 
following problems were encountered in using the SVDS LEA processor 

a. the LEA processor can read only one input tape 

b. the processor is limited to reading 500 files of data. 

The limitation of one input data tape for LEA seeins to imply that 
the nomiriil arid ;.ll e. - ir source trajectory sii'ulations must be 
run in ore SVDS jco ,> .c.. . . ■’!. The UfilVAC systc.Pi does not allow 


reading and then writing on the same tape, so It Is not possible to 
create a nominal trajectory tape In one job submittal and 
subsequently add data files for the erroi^ source cases as they are 
generated In other job submittals. However, Investigation of 
UNIVAC tape capabilities reveals that several existing tapes can be 
combined by executing a tape processor called MATCH. This allows for 
dividing the SVDS trajectory simulations Into several small jobs (each 
creating a tape of dispersed trajectory parameters) and combining 
the tapes for LEA purposes. The control card setup for this opera- 
tion is 

XQT MATCH 

DUP unltp unItQ 
TEF unltjK 

DUP un1t2, unlt^ 

TEF unitj. 


DUP unit^, unitQ 
TEF 

where unit^., i = 1, 2, . . . are the SVDS tapes with the dispersion 
data from the error source cases and unit^ is the tape unit onto 
which the files arc being combined. The TEF directive is required 
to separate the files of unit^ by an end of file (FOJ). 
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The process can be used to delete a file if a series of SVDS simula- 
tions were generated and one case was found to be in error. The 
MATCH directive will allow for submittal of several short SVDS simu- 
lations rather than one large multi-case simulation requiring a one 
hour run estimate. 

As previously indicated, when generating the SVDS trajectory simu- 
lations tape records should be written only at the flight events for 
which LEA processing is desired. This restriction is imposed by the 
limitation of storage location for the LEA processor. The processor is 
currently limited to 500 storage points of trajectory data. That 
is, the total number of the event and time slices for all perturbation 
cases and the nominal must not exceed 500: 

(2 Events + STiine slices)* 

(Perturbation cases + 1) £ 500. 

This limitation may not be realistic for dispersion analyses of 
trajectories such as reference mission 3A where an abort once around 
is considered. 

2.3.3 Conclusions of LEA Investigation 

As part of this study, the LEA computations were verified by hand 
calculating the following: 

a. the rotation matrix used to rotate ECI state vector deviations 
into the LI!S. 

b. the RSS of tfio slate vector deviations. 
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c. the covariance matrix at randomly selected points. 

The hand calculations match the computer determined data. 

The major difficulty encountered while attempting to execute the LEA 
processor was the lack of documentation defining input procedures 
and requirements. SD6 indicates that updated LEA processor docu> 
mentation is underway. 

3.0 CONCLUSIONS 

The SVDS program is now an adequate simulation tool for conducting 
error analyses of IMU errors. SVDS needs to be modified before a 
complete dispersion analysis can be conducted. In particular, SVDS 
atmospheric perturbation capability needs to be corrected, core 
requirements need to be reduced, and the problems of stacking cases 
which require tabular input data need to be resolved. The LEA 
processor is executing properly but better documentation is required 
to make it an effective tool. 
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IW ERROR SOURCE DEFINITION 





GAP/SVOS TPAJECTORY DATA COWARISOK 




CA«»/SyOS TRAJECTORY DATA COHPARISOA 




. table III (continued) 
6AP/SV0S TRAJECTORY DATA COMPARISON 
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